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ABSTRACT 


This  project  investigates  the  problem  of  locating  seismic  events  from  combined  data  sets  of  regional  and  teleseismic 
arrival  times,  based  on  the  use  of  a  unified  3D  model  of  the  Earth's  velocity  structure  to  predict  travel  times  for  both 
types  of  data.  Inherent  to  this  problem  is  the  joint  tomographic  calibration  of  the  Earth  model  with  both  regional  and 
teleseismic  ground-truth  data.  The  goal  of  the  unified  approach  is  to  remove  the  inconsistencies  that  result  when 
travel-time  predictions  are  performed  with  a  mixture  of  separately  calibrated  regional  and  global  models  and 
empirical  corrections,  which  can  lead  to  degradation  of  event  location  performance.  Numerical  experiments  we  have 
performed  with  joint  regional/teleseismic  travel-time  tomography  have  demonstrated,  for  example,  that  joint 
calibration  preempts  a  baseline  shift  between  regional  and  teleseismic  travel-time  predictions. 

Our  project  has  focused  largely  on  the  practical  difficulties  of  implementing  a  unified  location/calibration  capability. 
One  such  difficulty  is  the  computational  challenge  of  performing  tomographic  inversion  with  very  large  numbers  of 
data  and  model  parameters.  Available  data  bases,  for  example,  contain  many  millions  of  regional  and  teleseismic 
arrivals  while  3D  Earth  models  with  appropriate  resolution  contain  millions  of  grid  nodes.  Travel-time  computation 
in  3D  models  presents  another  significant  challenge.  It  can  require  days  of  central  processing  unit  (CPU)  time  to 
forward  model  a  travel-time  data  set  when  a  hilly  3D  raytracing  methods  is  used. 

An  efficient  alternative  to  full  3D  raytracing  is  travel-time  linearization,  which  approximates  the  travel  time  in  a  3D 
model  by  integrating  its  slowness  function  along  the  raypath  computed  in  a  ID  reference  model.  Our  numerical 
experiments  with  this  approach  have  shown  that  it  is  quite  accurate  for  teleseismic  travel  times,  as  judged  by 
comparing  to  travel  times  calculated  with  a  highly  accurate  raybending  method.  While  we  previously  had  dismissed 
linearization  as  being  inadequate  for  calculating  travel  times  at  regional  distances,  we  have  recently  pursued  a 
modification  of  the  linearization  technique  that  shows  promise  at  least  for  distances  up  to  ten  degrees.  The 
modification  involves  warping  of  the  Mohorovicic  (Moho)  discontinuity  and  other  velocity  interfaces  in  the  3D 
model  to  match  their  flat  counterparts  in  the  ID  reference  model.  This  paper  explains  the  warping  technique  and 
presents  some  examples  showing  its  promise  for  modeling  regional  travel  times.  We  also  present  preliminary  results 
of  joint  regional/teleseismic  tomography  based  on  the  new  linearization  method. 
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OBJECTIVE 


Several  studies  have  shown  that  3D  models  of  crust  and  upper  mantle  velocity  structure  predict  significantly  more 
accurate  travel  times  for  regional  seismic  phases  than  do  ID  global  models,  leading  to  more  accurate  locations  for 
small,  regionally  recorded  events  (e.g.  Ryaboy  et  al.,  2001;  Reiter  et  ah,  2005;  Murphy  et  ah,  2005;  Flanagan  et  ah, 
2007).  However,  it  is  not  uncommon  to  also  observe  teleseismic  arrivals  from  small  events  of  monitoring  interest,  at 
least  in  part  of  the  teleseismic  distance  range.  While  adding  teleseismic  to  regional  data  should  in  theory  improve 
event  locations  even  further,  it  has  not  always  been  the  case  in  practice.  For  example,  Yang  et  ah  (2004)  found  that 
event  locations  obtained  by  combining  regional  and  teleseismic  data  were  only  more  accurate  than  locations  found 
from  the  separate  data  sets  if  they  corrected  for  a  bias  between  teleseismic  and  regional  travel-time  predictions,  each 
made  from  a  different  3D  model:  CUB1.0  (Ritzwoller  et  ah,  2002)  for  regional,  and  J362D28  (Antolik  et  ah,  2003) 
for  teleseismic.  The  correction  they  inferred  is  small  (0.79  s)  but  even  small,  systematic  shifts  in  teleseismic  travel 
times  can  induce  large  location  errors  owing  to  the  high  apparent  velocity  of  teleseismic  arrivals. 

We  are  investigating  a  methodology  for  seismic  event  location  based  on  the  principle  that  travel  times  for  all  seismic 
phases  should  be  predicted  consistently  with  a  single  3D  Earth  model.  This  unified  modeling  approach  would  be 
applied  in  locating  events  from  combined  sets  of  regional  and  teleseismic  data,  and  in  a  joint  tomographic  inversion 
to  calibrate  the  3D  model.  The  feasibility  and  value  of  the  unified  3D  modeling  approach  has  been  demonstrated 
recently  by  Myers  et  ah  (2011)  and  Simmons  et  ah  (2011)  in  companion  papers.  Our  efforts  in  this  project  have 
focused  on  developing  computational  methods  that  improve  the  efficiency  and  practicality  of  joint 
regional/teleseismic  location  and  calibration. 

RESEARCH  ACCOMPLISHED 


Our  project  is  focusing  on  the  development  of  a  prototype  capability  for  event  location  and  travel-time  calibration 
with  3D  Earth  models,  applicable  to  regional  and  teleseismic  data.  Our  goal  is  not  to  develop  a  complete  and 
operations-ready  3D  capability,  but  rather  to  demonstrate  the  value  and  practicality  of  the  unified  approach.  The 
following  sections  describe  our  progress  to  date  in  developing  two  key  components  of  the  approach:  3D  travel-time 
calculation  and  joint  regional/teleseismic  tomography. 

Travel-Time  Linearization 


A  number  of  numerical  methods  are  available  for  raytracing  and  travel-time  calculation  in  3D  Earth  models,  such  as 
the  finite-difference  eikonal  method  (e.g.,  Podvin  and  Lecomte,  1991),  fast  marching  (Rawlinson  and  Sambridge, 
2004),  and  two-point  raybending  (e.g.,  Um  and  Thurber,  1987).  Given  the  computational  intensity  of  these  full  3D 
methods,  many  applications,  such  as  global  tomography  with  teleseismic  data,  have  used  approximate  travel  times 
calculated  with  linearization  around  a  ID  reference  model.  We  will  describe  this  approach  as  it  is  conventionally 
applied  and  then  propose  a  modification  of  the  approach  that  attempts  to  improve  its  accuracy. 


Let  r  denote  the  position  vector  in  some  coordinate  system  (e.g.,  with  origin  at  the  Earth's  center),  and  consider  the 
travel  time,  T,  between  a  source  point,  rj,  and  receiver  point,  r2,  through  a  3D  velocity  function,  v(r).  Geometrical 
ray  theory  implies  that  T  is  obtained  by  integrating  the  slowness  function,  l/v(r),  along  the  geometrical  ray 
connecting  rj  and  r2.  Denoting  the  geometric  ray  as  F,  we  can  express  this  as 


T  = 


•'ri  v(r)  1  ” 


(1) 


where  | dr | denotes  the  differential  distance  along  the  ray.  Fermat's  Principle  states  that  T  is  stationary  with  respect  to 
perturbations  to  the  ray  path,  allowing  us  to  approximate  the  travel  time  as 


T 


Jr24-|df|, 

Jri  r(r) 


(2) 


where  r  is  an  approximate  ray  connecting  rjand  r2.  The  accuracy  of  the  approximation  is  determined  by  the 
proximity  of  r  to  r. 
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The  standard  approach  to  travel-time  linearization  applies  Equation  2  with  r  taken  as  the  geometrical  ray  in  a 
reference  velocity  function,  vref.  Denoting  this  ray  as  rref,  Equation  2  becomes 


ri  v(rref) 


|rfrref  I- 


(3) 


The  approximation  becomes  advantageous  when  vref  is  chosen  such  that  the  computation  of  rref  is  more  efficient  than 
the  computation  off,  e.g.  when  vref  is  a  ID  velocity  model. 


Example 

In  the  previous  year  of  the  project  we  conducted  numerical  experiments  to  test  the  accuracy  of  travel-time 
linearization  for  computing  travel  times  in  an  early  version  of  Lawrence  Livermore  National  Laboratory’s  (LLNL) 
3D  Earth  model,  LLNL-G3D  (Simmons  et  al.,  2009).  The  experiments  considered  travel  times  between  four  station 
locations  and  a  3D  grid  of  event  locations.  The  event  grid  covered  the  box  15°-55°N,  30°-80°E  (roughly,  the 
Middle  East)  at  1°  spacing  and  contained  five  event  depths  between  the  surface  and  120  km.  The  four  stations,  each 
at  zero  depth,  were  located  at  the  corners  of  the  event  grid.  Travel  times  for  the  resulting  41,820  event-station  paths 
were  calculated  with  the  standard  linearization  method  using  the  AK135  ID  model  (Kennett  et  al.,  1995)  as  the 
reference  model.  These  times  were  compared  to  highly  accurate  raybending  calculations  through  LLNL-G3D 
performed  with  the  algorithm  described  by  Simmons  et  al.  (2009). 
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Figure  1.  Discrepancy  between  travel  times  calculated  through  model  LLNL-G3D  with  two  methods: 

raybending  and  linearization.  Linearized  times  were  calculated  with  AK135  as  the  reference  model. 
Each  panel  shows  the  discrepancies  for  paths  connecting  a  geographic  grid  of  events  at  15-km  depth 
to  one  of  four  surface  stations  at  the  corners  of  the  event  grid.  The  travel-time  scale  has  been  clipped 
at  ±0.7  s. 
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Some  results  of  these  tests  (also  shown  in  Rodi  et  al.,  2010)  are  shown  in  Figures  1  and  2.  The  figures  display  the 
discrepancy  between  linearized  and  raybending  times,  which  can  be  interpreted  as  the  error  in  the  linearization 
calculations  given  the  high  accuracy  of  the  raybending  method.  Figure  1  plots  the  travel-time  discrepancies  as  a 
function  of  event  epicenter,  for  the  events  at  15-km  depth  and  for  each  of  the  four  corner  stations  (labeled  NW,  NE, 
etc.).  Figure  2  displays  the  same  results  as  a  histogram  binned  by  travel-time  discrepancy  and  epicentral  distance. 
These  results,  and  those  for  the  other  event  depths,  indicate  that  the  errors  in  linearized  times  at  teleseismic  distances 
are  generally  well  below  0.5  s.  At  regional  distances,  however,  the  errors  can  be  several  seconds  or  more  (up  to  30  s 
for  some  paths).  A  major  contributor  to  the  large  errors  at  regional  distances  is  the  wide  variability  of  Moho  depth  in 
Eurasia,  from  about  10  km  under  the  Indian  Ocean  to  75  km  beneath  the  Himalayan  Mountains.  Thus,  for  example, 
some  first-arriving  Pn  rays  in  LLNL-G3D  may  be  contained  in  the  AK135  crust. 
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Figure  2.  Discrepancy  between  travel  times  calculated  through  model  LLNL-G3D  with  two  methods: 

raybending  and  linearization.  The  travel-time  discrepancies  shown  are  the  same  as  in  Figure  1,  but 
displayed  as  histograms  binned  with  respect  to  time  discrepancy  (with  0.2  s  increment)  and 
epicentral  distance  (1°  increment).  The  time  axis  is  clipped  at  7.5  s. 


Finearization  with  Ray  Warping 

The  modification  of  the  linearization  method  we  are  investigating  attempts  to  find  a  better  approximation  to  the 
geometrical  ray,  F,  in  the  3D  velocity  model,  v(r).  It  assumes  the  approximate  ray  has  the  form 
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f  =  W  (r'ref),  (4) 

where  W  is  a  specified  warping  transformation.  The  ray  r'ref,  like  rref,  is  calculated  in  vref  but  its  endpoints  are 
modified  by  the  warping  transformation  to  be 

r'  =  W~\  r,)  (5) 

r2  =  W  1(r2),  (6) 

where  W_1  is  the  inverse  transformation  of  W.  This  ensures  that  r  has  the  correct  endpoints,  rx  and  r2.  Substituting 
Equation  4  into  Equation  2  yields 


rr2  1 

Jri  v(W( fr'ef)) 


\dW  (rr'ef)|. 


(7) 


This  integration  along  the  warped  reference  ray  can  be  converted  to  integration  along  the  reference  ray  itself  with  a 
change  of  variable.  Define  the  quantity  y  at  each  point  of  a  geometrical  ray  r  as 


7(f)  = 


|dW(f)| 

|df| 


(8) 


The  quantity  y  is  the  ratio  between  the  differential  distance  along  the  ray  before  and  after  warping;  it  is  determined 
by  the  local  direction  of  the  ray  and  the  Jacobian  of  W.  The  conversion  of  Equation  7  to  an  integral  along  the 
unwarped  reference  ray  then  becomes 


Q 


r2  r(rref) 


i  v(W( rref)) 


|dr, 


ref  I 


(9) 


The  objective  in  this  approach  is  to  choose  W  so  that  it  transforms  the  reference  ray  into  a  better  approximation  of 
the  geometrical  ray  in  the  3D  model.  This  is  consistent  with,  if  not  equivalent  to,  transforming  the  reference  model 
into  a  better  approximation  of  the  3D  model,  i.e.,  making  vref  (W-1(r))  resemble  v(r)  more  closely  than  vref(r)  does. 
Our  initial  experiments  with  the  approach  have  attempted  to  achieve  this  by  choosing  W  to  move  points  vertically  in 
such  a  way  that  the  velocity  discontinuities  of  the  reference  model  and  3D  model  coincide.  Then,  for  example, 

W  (rr'ef)  will  be  at  the  Moho  discontinuity  of  the  3D  model  when  rr'ef  is  at  the  Moho  of  the  reference  model,  ensuring 
that  a  Pn  ray  in  the  reference  model  traverses  the  upper  mantle  of  the  3D  model. 


Example 

To  test  whether  ray  warping  improves  the  accuracy  of  linearized  travel-time  calculations,  we  repeated  the  numerical 
tests  with  model  LLNL-G3D,  described  above,  replacing  the  standard  linearization  method  with  the  new  approach. 
We  defined  the  warping  transformation  so  that  the  depths  of  all  interfaces  in  the  ID  reference  model  matched  those 
of  LLNL-G3D,  which  include  the  ocean  bottom;  interfaces  between  as  many  as  three  sedimentary  layers  and  three 
deeper  hard-rock  crustal  layers;  and  two  mantle  discontinuities  (at  410  and  660  km  in  AK135).  This  required  adding 
thin  (1-meter)  water  and  sedimentary  layers  to  the  AK135  profile  to  make  the  reference  model  compatible  with 
LLNL-G3D. 

Figure  3  shows  the  travel-time  discrepancies  between  linearized  travel-time  calculations,  performed  with  ray 
warping,  and  raybending  calculations.  The  figure  displays  the  discrepancies  as  distance-dependent  histograms  and  is 
analogous  to  Figure  2.  Comparing  Figures  2  and  3,  we  can  see  that  ray  warping  greatly  reduces  the  errors  in 
linearized  times  at  regional  distances,  up  to  about  10  degrees  of  so.  The  errors  at  teleseismic  distance,  however,  are 
increased  slightly.  Ray  warping  does  not  reduce  the  errors  at  epicentral  distances  between  about  14°  and  24°  to 
acceptable  levels.  This  distance  range  corresponds  to  AK135  rays  that  turn  in  or  near  the  transition  zone. 

Figure  4  shows  the  results  of  additional  accuracy  tests  involving  events  at  four  other  depths.  The  discrepancies  in 
these  cases  are  also  much  smaller  at  regional  distances  than  with  standard  linearization,  although  some  problems 
occur  for  surface  events  (which  include  events  at  the  top  of  the  ocean),  and  the  deepest  events,  whose  rays  are 
affected  by  the  transition  zone  at  shorter  distances. 
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Figure  3.  Discrepancy  between  travel  times  calculated  through  model  LLNL-G3D  with  two  methods: 

raybending  and  linearization  with  ray  warping.  Linearized  times  were  calculated  with  AK135  as  the 
reference  model.  The  discrepancies  are  displayed  as  histograms  binned  by  time  discrepancy  and 
epicentral  distance. 


Joint  Regional/Teleseismic  Tomography 

Rodi  et  al.  (2010)  showed  preliminary  results  of  applying  joint  regional/teleseismic  tomography  to  update  the  crust 
and  upper  mantle  P-wave  velocity  model  for  southern  Asia  developed  by  Reiter  et  al.  (2009).  That  model,  called 
JWM,  was  derived  with  the  joint  inversion  of  regional  (predominantly  Pn)  travel  times,  which  determined  the  P 
velocity  structure,  and  Rayleigh-wave  group  delays,  which  determined  the  S  velocity  structure,  applying 
geostatistical  constraints  to  couple  the  two  data  sets.  We  updated  the  JWM  P-wave  model  by  applying  joint 
tomography  to  the  same  Pn  data  set  used  to  derive  JWM  combined  with  teleseismic  travel  times.  Only  a  single  step 
of  an  iterative  inversion  was  done,  i.e.  linearized  tomography  was  performed.  The  resulting  model  fit  the  regional 
travel-time  data  as  well  as  JWM  while  also  removing  a  bias  of  over  one  second  in  the  teleseismic  times  predicted  by 
JWM,  thus  demonstrating  the  merits  of  joint  calibration. 

We  are  redoing  our  joint  regional/teleseismic  tomography  exercise  to  reflect  two  enhancements  we  made  to  our 
tomography  algorithm.  First,  we  expanded  the  capacity  of  the  algorithm  to  handle  up  to  one  million  travel-time  data 
and  one  million  velocity  model  parameters  (running  on  a  serial  processor  accessing  16  Gbytes  of  memory),  which 
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allowed  us  to  include  data  from  both  teleseismic  events  and  teleseismic  stations,  together  with  the  regional  times,  in 
one  inversion.  Previously  we  found  it  necessary  to  treat  the  two  types  of  teleseismic  paths  separately.  The  final  data 
set  comprises  about  140,000  regional  arrivals  (those  used  to  derive  JWM)  and  660,000  teleseismic  P  arrivals 
covering  event-station  distances  from  25°  to  55°. 
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Figure  4.  Discrepancy  between  travel  times  calculated  through  model  LLNL-G3D  with  raybending  and 

linearization  with  ray  warping,  displayed  as  histograms  binned  by  time  discrepancy  and  epicentral 
distance.  Linearized  times  were  calculated  with  AK135  as  the  reference  model.  Each  panel  displays 
the  discrepancies  for  paths  from  Station  SE  to  a  geographic  grid  of  events  at  the  labeled  focal 
depth,  z. 


The  second  enhancement  we  made  to  our  tomography  algorithm  was  the  incorporation  of  our  new  method  for  travel¬ 
time  linearization,  based  on  ray  warping,  as  an  optional  forward  model  for  regional  and  teleseismic  data.  Previously, 
the  algorithm  assumed  that  regional  times  were  linearized  (in  the  standard  way)  with  respect  to  a  3D  model  and  that 
teleseismic  times  were  linearized  relative  to  AK135.  Normally,  the  3D  reference  model  for  regional  times  is  the 
model  resulting  from  the  previous  iteration  of  a  nonlinear  inversion  but,  in  this  exercise,  it  is  JWM.  We  use  the 
Podvin-Lecomte  finite  difference  method  (Podvin  and  Lecomte,  1991)  to  compute  regional  travel  times  in  a  3D 
reference  model,  and  their  sensitivities  to  velocity  parameters.  Our  second  enhancement  allows  regional  times  to  be 
linearized  with  respect  to  AK135,  as  is  done  with  teleseismic  travel  times.  Additionally,  linearization  around  AK135 
can  be  done  with  or  without  ray  warping. 


As  a  prelude  to  generating  a  new  tomographic  update  of  JWM,  we  tested  the  accuracy  of  travel-time  linearization, 
relative  to  AK135,  for  modeling  the  regional  data  set  from  southern  Asia.  The  tests  compared  the  linearized  times 
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through  JWM  to  Podvin-Lecomte  calculations.  While  the  Podvin-Lecomte  method  is  not  as  accurate  as  raybending, 
it  is  still  a  full  3D  method.  Figure  5  shows  the  result  of  two  tests.  In  each,  linearized  times,  obtained  using  AK135  as 
the  reference  model,  are  compared  to  the  Podvin-Lecomte  times.  The  discrepancies  shown  in  the  left  panel  were 
obtained  with  the  standard  linearization  method,  while  those  on  the  right  used  linearization  with  ray  warping.  As  in 
the  earlier  example  with  LLNL-G3D,  we  see  that  the  use  of  ray  warping  dramatically  reduces  the  errors  in  regional 
travel-time  calculations.  The  errors  are  somewhat  larger  than  before  (c.f.  Figure  3),  but  errors  in  the  Podvin-Lecomte 
times  may  be  partly  to  blame. 
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Figure  5.  Discrepancy  between  travel  times  calculated  through  model  JWM  with  Podvin-Lecomte  raytracing 
and  linearization  around  AK135.  The  discrepancies  are  shown  as  histograms  binned  by  time  and 
epicentral  distance  (using  bin  sizes  of  0.1  s  and  1°,  respectively).  Left:  The  standard  linearization 
method  was  used.  Right:  Linearization  with  ray  warping  was  used. 

Figure  6  shows  a  vertical  slice  through  models  obtained  with  linearized  joint  inversion  of  the  regional  and 
teleseismic  travel-time  data  sets  from  southern  Asia.  For  both  models,  teleseismic  times  were  linearized  with  respect 
to  AK135.  Ray  warping  was  used,  although  for  teleseismic  distances  it  does  not  have  a  big  effect.  For  the  inversion 
model  in  the  top  panel,  regional  times  were  linearized  with  respect  to  the  3D  JWM  model,  as  had  been  done  in  our 
previous  experiments  with  joint  regional/teleseismic  tomography.  In  this  case,  Podvin-Lecomte  finite-difference 
raytracing  was  used  to  obtain  the  reference  times  through  JWM  and  their  sensitivities  to  velocity  parameters.  For  the 
inversion  model  in  the  bottom  panel,  regional  times  were  linearized  with  respect  to  AK135,  using  the  ray  warping 
technique  for  improved  accuracy.  We  see  that  the  models  are  very  similar  to  one  another,  further  supporting  the 
adequacy  of  linearization  of  regional  travel  times  with  respect  to  a  ID  reference  model.  Further,  the  main  features  of 
both  models  are  similar  to  those  of  our  earlier  joint  inversion  models  (Rodi  et  al.,  2010),  which  were  derived  from 
subsets  of  the  teleseismic  data. 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  numerical  techniques  we  have  developed  to  date  have  the  potential  to  significantly  enhance  the  practicality  of 
joint  regional/teleseismic  location  and  calibration  with  a  unified  3D  Earth  model.  In  particular,  our  modification  of 
travel-time  linearization  with  the  use  of  ray  warping  seems  to  extend  the  applicability  of  linearization  as  an  efficient 
alternative  to  3D  raytracing  to  regional  distances.  For  which  types  of  3D  structures  and  which  epicentral  distances 
and  event  depths  it  can  be  applied,  with  acceptable  and  predictable  accuracy,  will  require  further  study  to  determine. 

Improving  the  efficiency  of  travel-time  calculation  is  important  for  both  event  location  and  travel-time  calibration. 
The  task  of  locating  events  on  a  routine  basis,  using  arrival-time  picks  from  a  large  base  of  stations  changing  over 
time,  is  greatly  simplified  if  travel  times  for  arbitrary  paths  can  be  calculated  quickly.  Computational  efficiency  is 
also  crucial  in  building  a  unified  3D  Earth  model  via  tomographic  inversion  of  large  travel-time  data  sets.  Among 
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other  reasons,  it  facilitates  the  challenging  task  of  computing  model  uncertainty.  Beyond  computational  speed,  the 
validity  of  a  linear  approximation  to  travel  times  is  valuable  in  and  of  itself,  allowing  iterative  inversion  schemes  to 
converge  faster  and  increasing  the  validity  of  uncertainty  measures  that  are  based  on  linear  inverse  theory. 
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Figure  6.  A  vertical  slice  (Line  A-B  in  map  to  the  right)  through  two  3D  P  velocity  models  obtained  with 

linearized  joint  inversion  of  regional  and  teleseismic  P-wave  travel  times.  Each  is  shown  in  terms  of 
its  velocity  deviation  from  model  JWM.  The  inversion  models  are  both  based  on  linearization  of 
teleseismic  travel  times  relative  to  the  AK135  model,  but  used  different  modeling  technique  for 
regional  (Pn)  times.  Top:  Regional  travel  times  are  linearized  with  respect  to  JWM,  with  reference 
travel  times  and  their  sensitivities  computed  with  the  Podvin-Lecomte  method.  Bottom:  Regional 
times  are  linearized  relative  to  AK135,  using  ID  raytracing  with  ray  warping  to  compute  reference 
times  and  sensitivities. 
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